Critical Fluctuation of Wind Reversals in Convective Turbulence 
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The irregular reversals of wind direction in convective turbulence are found to have fluctuating 
intervals that can be related to critical behavior. It is shown that the net magnetization of a 2D Ising 
lattice of finite size fluctuates in the same way. Detrended fluctuation analysis of the wind reversal 
time series results in a scaling behavior that agrees with that of the Ising problem. The properties 
found suggest that the wind reversal phenomenon exhibits signs of self-organized criticality. 



In turbulent thermal convection at high Rayleigh num- 
bers (Ra) it has recently been observed that there exists 
not only large-scale circulating motion, called mean wind, 
but also abrupt reversals of the wind direction, whose 
physical origin is still largely unknown 0-0- For our 
purposes, Ra is simply a non-dimensional measure of the 
temperature difference between the bottom and the top 
plates of the container within which the convective mo- 
tion occurs. Metastable states have been suggested to 
describe the two opposite directions of the wind, and the 
reversal of its direction is to be understood in terms of 
the imbalance between buoyancy effects and friction 0- 
Instead of searching for the origin of the wind reversals 
in the framework of hydrodynamical considerations, we 
investigate in this paper the possibility of understanding 
the phenomenon in a totally different context, namely: 
critical phenomenon. We shall find a measure to quantify 
the fluctuations in the wind direction, and then demon- 
strate that its behavior corresponds to one exhibited by a 
system undergoing a second-order phase transition. We 
then perform a detrended fluctuation analysis to deter- 
mine the detailed properties of the fluctuations of the 
wind, more specifically its scaling behavior. 

The experimental data that we analyze are the same 
as those reported in |^] and studied in 7]. By varying 
the pressure and lowering the temperature of the gas, 
the Rayleigh number could be varied between 10^ and 
10^^. Further details of the apparatus can be found in 
0- We focus on the data that give the wind speed and 
direction for a continuous period of up to one week at 
Ra= 1.5 X 10^^. Figure 1 shows a small segment of the 
wind velocity data for 6.5 hr, starting at an arbitrary 
time. Note how the wind changes direction suddenly in 
the time scale of that figure. We proceed directly to an 
interpretation of the fluctuations of the wind velocity. 

In a fully developed turbulent convection at high Ra 
there are two opposing dynamical features. One is the 
emission of plumes from the top and bottom boundary 
layers; they occur at random locations and at random 
times in varying sizes. The other is the existence of mean 
wind that rotates in one direction or another, making 
rapid reversals at seemingly random intervals. We regard 
the former as the disordered motion of the components 
of a complex system, and the latter as the ordered mo- 



tion of the whole of the system. For low Ra (say below 
10^), the ordered motion is not sufficiently impeded by 
the disordered motion to cause reversals of the wind di- 
rection. At high Ra, the cumulative effect of the many 
plumes that is strong enough to reverse the wind direc- 
tion. The system then proceeds as before except that the 
wind rotates in the opposite direction with varying mag- 
nitude until another reversal occurs due to the collective 
action of the disorganized plumes. If the system is at a 
critical state, whether self-organized or not, the competi- 
tion between the ordered and disordered motions leads to 
the wind switching directions at irregular intervals of all 
scale. The probability of occurrence of the wind duration 
T between reversals should satisfy a power law 



p{t) 



(1) 



as a manifestation of criticality. Such a power law has 
been found in the data 0- The discussion above de- 
scribes our view of the origin of such a scaling behavior. 

We now advance the idea that the above description of 
the wind and plumes in convective turbulence in terms 
of ordered and disordered motions has its corresponding 
counterparts in the 2D Ising model of critical behavior. 
In the Ising system of near-neighbor interactions without 
external magnetic field the lattice spins tend to align in 
the same direction except for the random disorientation 
due to thermal fluctuation. For a finite lattice the net 
magnetization, M, is non- vanishing. For T < Tc, the 
critical temperature, M is likely to persist in the same 
direction for longer time in lattice-spin updating than at 
higher T . At T > the thermal interaction dominates, 
and M is more likely to flip sign more frequently upon 
updating. The fluctuation of the signs of M is therefore 
a property that reflects the tension between the ordered 
and disordered interactions of the whole system. 

Since the mean wind is a global phenomenon in a vessel 
of finite volume, it is sensible for us to associate the wind 
direction with the sign of M of the Ising lattice of finite 
size. We can then map wind reversal to the reversal of 
M upon updating the lattice spins in a simulation. The 
plumes are the disordered fluctuations that correspond to 
the spin fluctuations due to thermal agitation, and the 
wind is the ordered motion that can change direction just 
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FIG. 1: A segment of the data on wind velocity fluctuation 
at Ra= 1.5 X lO" . The data were obtained from an apparatus 
that was essentially a cylindrical container of 50 cm diameter 
and 50 cm height, filled with cryogenic helium gas. Two sen- 
sors were placed outside the boundary layer on the sidewall 
of the container near the middle section of the container, and 
were separated vertically by a distance of 1.27 cm. 



as the magnetization can change sign when enough lat- 
tice spins change directions. The key connection between 
the two problems is the mapping of the real time in tur- 
bulence to the time of updating the Ising configurations. 
It is therefore crucial that each configuration has some 
memory of the previous configuration before updating; 
hence we employ the Metropolis algorithm, which does 
precisely this. It should be noted that we are entering 
into a rather unexplored territory where the process of 
computer simulation itself is endowed with some physical 
significance, quite unrelated to the large body of analyt- 
ical work that has been devoted to the Ising model of 
infinite lattice. Our task is to show that the wind re- 
versal problem corresponds to the Ising problem of finite 
lattice at the critical temperature. 

To be more specific, we consider a square lattice of size 
L^, where L is taken to be 255, an odd number. We start 
with the site spins having a random distribution of ±1 
values. We then visit each site and determine from the 
usual near-neighbor interaction whether its spin should 
be reversed: yes, if the energy is lowered by the flip; if 
not lowered, the flip can still take place according to a 
thermal distribution specifled by temperature T. One 
time step is taken by the whole system when all sites 
are updated. We take 3 x 10^ time steps in total, and 
divide the whole series into 30 segments. The values of 
M at each of the 10** time points in each segment are 
discretized to ±1, according to M^O. A continuous string 
of M of one sign, either -1-1 or —1, forms a duration that 
is analogous to the mean wind rotating in one direction. 
The reversals of M correspond to the reversals of wind. 
Near the critical point, durations of all lengths can occur. 



Before considering the issue of criticality for a finite lat- 
tice, let us discuss the measure that we shall use for quan- 
tifiying the duration fluctuations appropriate for both 
the wind and Ising problems. The experimental data 
on wind consist of 8 segments, each having T = 10, 282 
time points. For the Ising case we have 30 segments, each 
having T = ID'*, roughly the same as wind data. Let N 
denote the number of reversals in a segment. With the 
locations of the reversals denoted by ti, i = 1, ■ ■ •, iV, de- 
fine Ti — ti^i — ti to be the ith duration (or gap), where 
to and tj^^i are assigned to be the left and right ends of 
the segment, respectively. Now, define the moment ^3 

1=0 

where q is any positive integer. Clearly, Gq — 1 and Gi = 
1/(A'' + 1). Gq is a measure that quantifies the pattern 
of reversals in each segment. For large q,Gq is a small 
number, since t^/T is small. Its value can be dominated 
by a few large gaps, as when T < Tc, or may become 
significant from the sum over many small contributions 
due to many small gaps, as when T > Tc- For a measure 
of the fiuctuations of Gq from segment to segment, we 
define an entropy-likc quantity 10^, IL] 

= -(G,lnG,), (3) 

where (• ■ •) implies an average over all segments. For 
brevity we shall refer to the study of the time series in 
terms of Sq as the gap analysis. In Fig. 2 we show by 
filled circles the result of the gap analysis on the wind 
data at Ra= 1.5 x 10^^. It is evident that for g > 2 the 
points can be well fitted by a straight line, shown by the 
solid line, exhibiting an exponential behavior for Sq 

In Sq = -Xq + Xo , A = 0.264 . (4) 

For the Ising simulation we must first decide on the 
proper value of the critical temperature Tc for a finite 
lattice. For an infinite lattice its value has been deter- 
mined analytically to be 2.269 in units of J/ks, where 
J is the coupling strength of near-neighbor interaction 
and the Boltzmann constant ^2 ■ For a finite lattice 
the value of Tc should be higher. We have performed 
the simulation of our Ising system at three values of T, 
and determined the properties of M reversal. In Fig. 
2 we show the results of our calculated values of Sq at 
T = 2.305,2.310 and 2.315. Only the one at T = 2.310 
(lowered by a factor of 2 for clarity) shows a nearly linear 
dependence in the plot. The dashed line is a linear fit of 
the points in open circle, giving a slope of A = 0.261. 
At the two neighboring values of T, the q dependencies 
of logSq (shown by triangles and squares) are not linear, 
the values at high q being higher than at T = 2.310. 
The linear behavior at T = 2.310 is almost the same 
as in the wind reversal problem, as can be seen visually 
by the dash-dot line, which is a parallel transport of the 
solid line for comparison, but displaced slightly from the 
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dashed line to avoid overlap. We regard T = 2.31 as the 
critical temperature in our Ising system, since it has 
the unique property of being different from those of the 
neighboring T on both sides. When T < Tc, the gaps are 
longer and Gq is larger at large q (but still ^1) with 
the consequence that Sq is larger. When T > Tc, the 
gaps are shorter, but many gaps can contribute in the 
sum in Eq. Q, resulting in Gq still being larger at large 
q with the consequence that Sq is also larger. It is only 
at the critical point that gaps of all sizes can occur, re- 
sulting in Gq to be smaller and therefore Sq also smaller 
at large q. Thus the exponential decrease of Sq is a sig- 
nature of criticality. The value of Tc obtained here is in 
accord with the result of another calculation, in which 
the normalized factorial moments are found to exhibit 
scaling behavior at the critical point, different from the 
non-scaling behaviors at neighboring T • In that cal- 
culation the measure studied quantifies the fluctuation of 
the cluster sizes in an Ising system on a square lattice of 
size L = 288, for which Tr is found to be 2.315. 
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FIG. 2: Moments in the gap analysis of wind reversal (filled 
circles) and magnetization reversal in Ising lattice (open sym- 
bols) for difi^erent temperatures. The open circles are lowered 
by a factor of 2 to give space for clarity. The solid line is a 
linear fit of filled circles, and the dashed line is a linear fit of 
open circles. The dash-dot line is parallel to the solid line, 
placed near the dashed line for comparison. 

The normalizations of Sq for the wind and Ising prob- 
lems are not the same, since the average numbers N of 
reversals are different. However, the exponential behav- 
iors are remarkably identical. The q dependence of Sq is 
a quantitative measure of the fluctuation behavior of the 
reversals. The fact that the slope A is the same for both 
the wind and magnetization problems suggests strongly 
that the wind reversal in convective turbulence at high 
Ra is a critical phenomenon. Moreover, since we have not 
tuned any adjustable parameter in the wind problem to 
bring the system to the critical point, as we have done for 
the Ising system by varying T, we conclude that the wind 
reversal phenomenon is a manifestation of self-organized 



criticality (SOC) [ij]. 

We now search for a power-law behavior that charac- 
terizes changes in the wind direction. (For other such 
efforts, see Q). Our method is the dctrendcd fluctua- 
tion analysis (DFA), which has been found to reveal the 
scale-independent nature of time series in a variety of 
problems, ranging from heartbeat irregularity [Ts^ and 
EEG 13 to economics 17]. In that analysis we look for 
scaling behavior in the RMS deviation of the wind veloc- 
ity from local linear trends. Given the time series of the 
wind velocity V{t) over a total range of T^ax, we divide 
it into B equal bins of width k, discarding the remainder 
%na.x — Bk. Let Vb{t) denote the linear fit of V{t) in the 
6th bin. The variance of the deviation of V{t) from the 
local trend, Vb{t), in bins of size k is defined by 



B ^ k 

fc=i t=ti 



(5) 



where ti = 1 + (6 — l)fc and t2 — bk, measured in units 
of At = 5 sec, so that the values of t are dimensionless 
integers that count the time points in the data. The goal 
is to study the behavior of the RMS fluctuations F{k), 
as k is varied. If there is no characteristic scale in the 
problem, then F{k) should have a scaling behavior 



F{k) cx k° 



(6) 



This power law cannot be valid for arbitrarily large k 
because the series V{t) is bounded, so for very large k 
the linear trend is just the V{t) = line, and the RMS 
fluctuation F{k) must become independent oik. Thus we 
expect Ini^(fc) to saturate and deviate from © at some 
large k. We note parenthetically that we have applied 
DFA to the unintegrated time series V{t), which is a 
departure from the usual practice. 
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FIG. 3: Scaling behaviors of F(A;) in DFA of wind reversal 
(filled symbols) and magnetization reversal in Ising lattice at 
the critical temperature (open circles). Lines are linear fits. 

In Fig. 3 we show F(A:) in a log-log plot for four equal 
segments of the complete wind data in solid symbols. 
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The segment segl is for time running from to 116,435 
s, corresponding to T^^ax — 23, 287; other segments all 
have the same length. We have limited the maximum 
bin size to 2,580, so that even for the largest bin the fluc- 
tuations can be averaged over 9 bins. Evidently, there 
is a good scaling for each segment. The points for seg3 
and seg4 are shifted upwards by the quantities indicated 
in order to give clarity without overlap. Note that the 
segl data do not have the same magnitude of F{k) as 
the other segments; yet the scaling exponents are essen- 
tially the same. The deviation from the straight lines at 
the upper end is the saturation effect already discussed. 
There is another short region of scaling with a higher 
slope at low k. It is a consequence of fluctuations of the 
velocity within one direction of the wind, whose presence 
is evident in Fig. 1. Since the critical behavior identified 
here refers to wind reversals, and not to fluctuations of 
the wind velocity within one direction, we should ignore 
the lower short scaling region. 

In the scaling region to which we pay attention here, 
the slopes are a=1.20, 1.20, 1.21 and 1.22, for segl to 
seg4, respectively. The deviations among the segments 
are obviously small. The average value is 

a =1.21. (7) 

This large value of a implies a smoother landscape com- 
pared to the rough time series of white noise that is char- 
acterized by complete unpredictability Indeed, the 
fluctuations of the wind reversal time series has gaps of 
all sizes, the signature of critical behavior that is char- 
acterized by 1/f noise Q]. It is interesting to compare 
our result with the properties of the power spectral den- 
sity for the velocity found in Ref. [igj where a scaling 
behavior is shown to exist with a slope roughly -7/5 (not 
by fitting) in the region —3 < log / < —1.8. That range 



of frequency corresponds to 4.1 < ln(l//) < 6.9. If we 
identify the values of k in DFA to the time scale 1/f, then 
that range of ln(l//) corresponds to the range of Infc in 
Fig. 3, in which wc find the scaling behavior with the 
exponent a given in Eq. jT)). That value of a is not too 
different from 7/5. The scaling behavior found in DFA 
uses shorter segments of the whole data and exhibits the 
power law more precisely, from which the value of a can 
be more accurately determined. 

We now apply DFA to the Ising problem. We consider 
10 segments of the M reversal time series of the Ising 
lattice set at Tc, each segment having lO'' time points. 
From the F{k) determined in each segment, we average 
over all segments and show the resultant dependence on 
k in Fig. 3 by the open circles. Clearly, the points can 
be well fitted by a straight line. The slope is 

aM = 1.22, (8) 

which is essentially the same as that in Eq. Q for wind 
reversal. With the equivalence of these two scaling be- 
haviors established, we have found stronger evidence that 
the wind reversal problem is a critical phenomenon. 

To summarize, we have studied the time series of wind 
reversal in convective turbulence by two methods (gap 
analysis and detrended fluctuation analysis) and applied 
the same methods to the time series of the reversal of the 
net magnetization of a 2D Ising lattice. The same results 
are obtained for both problems. We therefore can assert 
that wind reversal exhibits all the essential properties 
characteristic of a critical behavior; apparently requiring 
no tuning, it can be regarded as self-organized. 
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